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Abstract 

In this paper we show that in addition to the known minimal surfaces which appear in 
the literature for computing the holographic entanglement entropy, there are other minimal 
surfaces with non-zero extrinsic curvature. We use the approach of regularization procedure 
presented by Fursaev et al in [7], to compute the quadratic and cubic curvature invariants 
on manifolds with squashed cones. The results can be used to find the leading and universal 
terms of the holographic entanglement entropy to understand which solution corresponds to 
the actual minimal entropy. 
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1 Introduction 

The entanglement entropy of a sub-system in a quantum field theory, is a non-local quantity 
which measures how this sub-system is correlated with its complement. A holographic de¬ 
scription of entanglement entropy first proposed in [lj and [2]. They argue that the leading 
divergence of the entanglement entropy of a region A in a CFT^+i is proportional to the area of 
a d-dimensional minimal surface in AdSd +2 whose boundary is given by dA. Generally one can 
find a series solution in terms of t/e for entanglement entropy, where t is the typical length of 
dA and e is the UV cut-off. In odd dimension d, there is a logarithmic divergence lnf/e, known 
as the universal term, which its coefficient is proportional to the central charges of dual CFT. 

Usually in computing the entanglement entropy in a quantum field theory one uses the replica 
trick by introducing a conical singularity. In the context of AdS-CFT one may use a similar trick 
by deforming the bulk metric so that the minimal surface contains a conical singularity. One of 
the recent works to derive the holographic prescription of [I] is the proposal of |3j. They show 
that in Einstein gravity the minimal surfaces can be found from the bulk equations of motion 
when we demand no singularity. 

To extend this prescription for higher derivative gravities several attempts have been done. 
For example in [3] it has been shown that the method used in [3] only works when the extrinsic 
curvature is small in Gauss-Bonnet gravity and results coincide with the known Jacobson-Myers 
entropy [5], which has established already as a correct entanglement entropy formula [6]. 

The other approach presented in [7], proposes a regularization procedure which has been 
developed to compute the integral curvature invariants on manifolds with squashed cones. By 
using this method one finds a holographic formula for entanglement entropy when the bulk grav¬ 
ity includes quadratic curvature terms. In special case of Gauss-Bonnet gravity this approach 
again gives the correct Jacobson-Myers entropy result. For further developments of this method 
to other theories of gravity see for example P ESI EH E2J E3 EH EH ESI EH EB] ■ 
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In this paper we have used the approach of |7] to find a formula for entanglement entropy of 
the six derivative gravity theories. Using our results and the known formula for the four deriva¬ 
tive gravity theories we will be able to study the holographic entanglement entropy of cylindrical 
and spherical surfaces which are living on the dual conformal held theories. We also compute 
the leading and the universal terms corresponding to the extremal surfaces with/without extrin¬ 
sic curvature. We will find two different types of solutions for these extremal surfaces, which 
exist in different regions of the parameter space of each theory. These regions can be found 
by demanding of positivity of the leading order terms in the entropy and reality condition of 
solutions. We will argue about how we can choose which solution gives the correct result i.e. 
the minimal entanglement entropy. 

Since we have used and developed the method of regularization procedure of Fursaev et al 
in [7] so let us review very briefly the main results of it: 

1. Consider a geometry M. and construct an orbifold M n made by cutting M along a 
co-dimension one hyper-surface in Ai, then glue n identical copies (replica method). The new 
geometry Ai n has a conical singularity located at r = 0 with the following metric (according to 
the replica method we need to set a period of 2tt n for angular coordinate r) 

ds 2 = r 2 dr 2 + dr 2 + r; x)dx l dx J , (1.1) 


where x l, s are the coordinates of the co-dimension two surface £ and 'jij is its intrinsic metric. 
This surface is chosen according to the entangling surface which we are interested to study in 
the dual conformal field theory. It can have different topologies in different dimensions. In 
this paper we have considered the cylinder and n-Sphere geometries embedded inside the AdS 
space-time. 

2. The geometry (11.11) includes the curvature singularity as well as the conical one. So in [7J 
a regularized metric (A4 n ) has been replaced by the above geometry as 

ds 2 = r 2 dr 2 + ^ ^ dr 2 + (a + r n c 1_n cos r) 2 ds|. (1.2) 

r 2 + b 2 

3. Inserting the regularized geometry into the integral curvature invariants gives rise to a 
result with an asymptotic series expansion in terms of parameter b when fe —>■ 0. Sending —>■ 1 
the result of each integral may contain a term which is proportional to 6°(1 — n) (i.e. independent 
of the regularization parameter). This limiting procedure in (7] suggests that one may perform 
the following replacement 

f dTdrd d ~ 2 Xy/~gR —>• n I drdrd d ~ 2 x\/QTZ + 47t(1 — n) f d d ~ 2 x^Ay + 0((1 — n) 2 ). (1.3) 
JM„ JM J X 


The left hand side describes the value of a curvature invariant term (here the Einstein-Hilbert 
term or in general any action constructed out of curvature invariants) on a regularized geometry 
A4 n , i.e. I[(AA n )\. The regularized geometry is given in equation (11.21) with periodicity r ~ 
r + 27rn. The right hand side contains two main terms. The first term I[Ai\ is equal to the 
same action as we begin with but computed on the non-regularized manifold A4 with metric 
given by (EH) and again r ~ r + 2irn. The second term is a contribution coming from the near 
region of the tip of the cone in the singular manifold M. This last term in above equation is 
noting but the area of the co-dimension two surface times the deficit angle of the cone. Now 
the holographic entanglement entropy can be computed by Shee = nd n (I[{M .„)] — I[M]) n ^i. 
This explains us why we need to expand I[(AA n )\ up to 0(n— 1). Explicitly starting from (11.31) . 
the holographic entanglement entropy will be proportional to the area of the co-dimension two 
surface which ends on the entangling region as conjectured by Ryu and Takayanagi. The same 
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computation has been done in [7] for quadratic curvature terms 


f drdrd d 2 Xy/gR 2 —>• n f drdrd d 2 xVQ1Z 2 

JMn J m 

+ 87 t (1 — to ) J d d ^ 2 Xy/^TZ + • • • , 

I dTdrd d ~ 2 XyfgR 2 —tn drdrd d ~ 2 xVQTZ 2 

Mn JM 

+ 4tt(1 - to) f d d ~ 2 x^f(n tw nfr^ - ^I< 2 ) + - 

J s 

dTdrd d - 2 Xy/~gRf iUa p —> n / dTdrd d ~ 2 xVQ'R 2 l 

m 4 M 


>2 


+ 8tt(1 — n) fd d 2 Xyfy(Uftvapnfnjnfn? - K^IC^) + 
J s 


(1.4) 


where n^’s are the unit mutually orthogonal normal vectors to £ (* = 1,2 for a co-dimension 
two surface) and K^ u are components of extrinsic curvature tensor defined as 

= K = (1.5) 


4. To find the holographic entanglement entropy S(E) for a surface £ in a quadratic curvature 
bulk gravity 

I[M] = — J d d Xy/g(^— + 2A + aR' 2 + bR 2 w + cR 2 „ a f^ , (1-6) 

we first compute the partition function. In semiclassical approximation the partition function 
is related to the action by — InZ = I[M\. In replica method we need to compute the value of 
7[A4 n ] by prescription mentioned above and then use the following relation 

Shee {£) = ( nd n — l)I[M n ] n -+\ . (1.7) 


Here the conical singularity of A4 n is located on a minimal co-dimension two hypersurface £ 
which ends on entangling surface £ on the boundary. Finally one gets the following formula for 
entanglement entropy 

S H ee {£) = 471- J d d ~ 2 x^i Q + 2aR + blZ^/n 1 - n 1 - + 2cT7 AtI/Q , (3 77,f n^n^vP- - ^K 2 - 2cK i flu Kp u ^j , 

( 1 . 8 ) 

where the four first terms are nothing but the Wald’s entropy and the last two terms are 
corrections due to existence of a non-vanishing extrinsic curvature of the minimal hyper-surface. 

The organization of the paper is as follows: In section 2 we review the holographic entangle¬ 
ment entropy in new massive gravity and find the differential equation which gives the minimal 
surface. We find two solutions for this equation with/without extrinsic curvature and compute 
the related leading terms and universal terms. In section 3 we add extensions to new massive 
gravity to study the behavior of the entanglement entropy in the presence of the new curvature 
cubed terms. To compute the integrals of curvature invariants we use the method reviewed 
above. In section 4 we revisit the general curvature square terms and in section 5 we do an 
exercise with quasi-topological gravity in five dimensions. In each section we discuss about the 
solutions in terms of their entanglement entropy in different regions of the parameter space. In 
last section we summarize and discuss about the results. 
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2 New Massive Gravity (revisited) 


Ini by using the method of squashed cones, the entanglement entropy in New Massive Gravity 
(NMG) [19J has been studied (see also m- They show that by using a bulk metric constructed 
using the Fefferman-Graham expansion, there exist a minimal surface and its corresponding 
entanglement entropy can be computed just by using the Wald’s entropy. Here in this section 
we show that in addition to the known minimal surface, there is another minimal surface with 
non-zero extrinsic curvature and we use the squashed cone method to compute the entanglement 
entropy. 

In three dimensions we consider a line as an entangling region, which lies on the boundaries 
of AdS 3 space-time located at z = 0. We Start from the AdS^ metric 

L 2 

ds 2 = g flu dx fJ ’dx l/ = — ( dt 2 + dr 2 + dz 2 ) , ( 2 . 1 ) 

where L is the AdS% radius. The Euclidean version of the NMG action is written as [2D] 


S 


wj ^ 3 


R + T 4AL 2 


R, IV R ,W - -R‘ 


( 2 . 2 ) 


where L = LyffoZ and satisfies 1 — + /£>A = 0, using the fact that (12.11) is a solution for 

the NMG equations of motion. To compute the entanglement entropy we use (USD, where for 
NMG it reduces to 


qNMG 

‘-’EE 




1 + 4A L 2 


TZ^n^n] 


-I*-* 


(2.3) 


To find the geometry of the extremal surface we must take t = 0 and r = f(z), then extremize 
(12.31) to find a suitable value for f(z). The strategy we will use through our computations is 
finding the tangent vectors at first and then making orthogonal unit vectors. In three dimension 
there is only one tangent vector which is given by e“ = dx a /dz = (0,1 ,f'(z)). Then we can 
build the induced metric of the minimal surface, 7 ^, and also its orthonormal normal vectors 


ds 2 = 7 ^dx^dx 1 * = 




ni^ = 


0 ,— 


Lf(z) 


L 


zVf'( z Y + 1 z Vf ( z ) 2 + 1 


n2 t i = 


L 

—, 0,0 

z 


(2.4) 


For each normal vector we have a corresponding extrinsic curvature. For example its 2:2: compo¬ 
nent is given by (the other non-zero components are proportional to this component when using 
equation (11.51) ') 


KL = 


L{f{zf-zf"{z) + nz)) 
z 2 (f'{z) 2 + 1 )§ 


KL = 0 . 


(2.5) 


The vanishing of K 2 Z is a consequence of having a Killing vector in the time direction. This 
is happening in all computations in this paper so we will drop the index i for simplicity in the 
up-comping computations. By inserting the above results into (12.31) . and after extremizing we 
find the following non-linear differential equation 

/ ooA (( 10/" 3 - 60 f' 2 f" 3 + 40 /" 7 ' 3 /" + 40 / ,/ 77 " - 4/^/ ,4 - 4 /( 4) - 8/ (4) /'V 

+ (30/ ,3 /" 2 + 30/7" 2 - 8 /'" - 8/"7 ' 4 - 16f"'f ,2 )z 2 + (2 - 4 /' 6 - 6/ ,4 )/'7 

- 2 f n - 6 f' 5 - 6/' 3 - 2 /') - (/' 6 + 3 /' 4 + 3 f 2 + 1 ) f"z + f 9 + 4 f 7 + 6/' 5 + 4 f 3 + f = 0 . 

( 2 . 6 ) 
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By a little investigation one can find two possible solutions for this differential equation 


fi(z) = \J z l - z 2 ; f 2 (z) = ^ z$ + 2z 0 qz - z 2 ] q 2 = 2\f ao - 1 , (2.7) 

where in both cases zq = f(z = 0) is equal to half of the length of the entangling line and 
the turning points are located at zt = zo and Zt = zo(q + y/\ + q 2 ) respectively (see figured]). 
Using the relation 1 — + f^X = 0 we see that to have a real valued solution we must restrict 

ourselves to either < 0 or > 2. In special value of foo = 2 both solutions coincide. The 
first solution has been found already in [9]. 



Figure 1: The dashed curve corresponds to fi(z) and the solid one is f 2 (z) with non-zero extrinsic 
curvature. 

The main difference between these two solutions is in the value of their extrinsic curvatures 

_ Lq(zQ + 2 z 0 qz — z 2 ) 


Kz 


/i( 2 ) 


= 0 , 


K 


h{z) 


z 2 zl(q 2 + 1)2 


( 2 . 8 ) 


Inserting the above results into (12.31) we can compute the entanglement entropy for both cases 


S { eI = S EE 


See = s ee 


/iO) 


ff 

^p J e 

27r L 

P-p V foo 


2vr 


dz 


Lz 0 


V( z o ~ z2 ) z2 


( 2 A/oo + 1 ) 


( 2 A/oo + 1 ) ( ln(—) + —— z 


4^2 
^ z 0 


zo 


(2.9) 


rzo 


mu t 


= —— I dz 


2LzqV2X 


p Je 


z V z o + 2qzz 0 — z 2 


47T 


= ZZlV 2A (ln(*») - -1 




20 


2 0 


+ 


( 2 . 10 ) 


Since the extrinsic curvature vanishes for fi(z), its corresponding entanglement entropy (|2.3D 
will be equal to the Wald’s entropy. As we can see for both solutions we have universal terms 

( 1 ) C i,._/A) x ci 2 tt L 27tL 3/qo — 2 


Q'- 1 ) = hi Ini —) li = 

EE 3 ( e j ’ 3 £ p vTo^ 


( 2 A/oo + 1 ) — 


3 

fSo 


S<?> = fa ln(—) . fa = ^LV2 a = 2 * L l° f °° - 1 


f 
J c 


2 

00 


( 2 . 11 ) 
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where c\ is the central charge of the dual CFT. For > 2 we have a real-valued function for 
f 2 (z) and both central charges are monotonically decreasing functions of /qq. In this region we 
have > c\ > C 2 > 0. Therefore the solution which minimizes the entropy is ,f 2 {z)- For 

2 > foo only the first solution exists as a minimal surface. In latter case demanding a unitary 
CFT dual restricts us to 2 > > |, see figure [2j 

Ci>0 c 2 >0 

! ■ , 

0 1 2 3 4 5 

/» 

Figure 2: The domain of validity (positive value of leading term in entropy) for different extremal 
surfaces in NMG. The blue domain, 2 > f 00 > §, corresponds to fi(z) and the red one, f aQ > 2, 
is for f 2 (z). 


3 Extended NMG 


To study the properties of entanglement entropy in presence of higher curvature corrections, in 
this section we consider the Extended NMG (ENMG), first introduced in [2D]. This is a theory 
which adds cubic curvature terms into the NMG and the action is given by 


S 



R +j2+ 4AL 



+ 1 V 

12 ^ 



72 rr + 64 r 
- rr 2 + - r 3 


(3.1) 


where for simplicity in notation we have defined R 2 = R ^ u RfiU and R 3 = R ^ a Ra ^ R ^. Similar 
to NMG the AdS 3 metric (12.111 is a solution of equations of motion for ENMG too. For ENMG 
we have 1 - /«> + /£>A + /£,// = 0. 

To compute the corrections to the entanglement entropy we need to find the corrections to 
the equation (12.3|) by using the method of squashed cones. In the presence of the 0(2) symmetry 
i.e. when the extrinsic curvature vanishes, we have the following relations for regularization of 
curvatures [8] 

R ^af} = + 2 7r ( 1 _ n )(n?n?Tqn? - nfnfn^n“)<5 s , 

Rf* = + 2vr(l - n)«<5 s , 

R = 1Z + 47t( 1 — n)<5s , (3-2) 


where 1Z curvatures are computed in the regular points A4 n /T, and 6% is the Dirac delta function. 
By combining these formulas and using the fact that ^ d d x ■ ■ ■ = n f M d d x • • •, we find 






d d Xy/g R 3 = 7i [ d d xVQTZ 3 + 67t(1 — n) [ d d 2 Xy/AfR^ a K va n^n IJ i +0((n- l) 2 ) , 

„ JM 4 s 

d d Xy/g RZ = n f d d xVQlZ i + 12ir(l — n) f d d ~ 2 Xy/^TZ 2 + 0((n — l) 2 ) , (3.3) 

f d d Xy/g RR 2 = n I d d xVGTZTZ 2 + 47t(1 - n) [ d d ~ 2 Xy/^ (JZ 2 + KR^n^n'') + 0((n - l) 2 ) , 
Jm„ jm 4s 

In absence of 0(2) symmetry, i.e. when the extrinsic curvature is not zero we need to consider 
the contribution of extrinsic curvature to the above relations. To do this we use the same steps 
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as JJj. We suppose that the possible independent allowed terms, according to the dimensional 
analysis, are those which are constructed either purely from extrinsic curvature or a combination 
of bulk curvature and extrinsic curvature. These are 

A' 4 , A'f, K 2 K 2 , KK 3 , K 4 . 

n^K^K, n^RPRi ;;, M 2 , M 2 , K ia , po RWR vo , (3.4) 

where 1Z,1Z^ U and IZ^ua /3 are defined on minimal co-dimension two hyper-surface £. For sim¬ 
plicity in use we call the first row as A"-terms and the second row AA"-terms. In the case of 
quadratic curvature terms, the cylindrical and spherical geometries in d = 4 were enough to 
compute the coefficients of R 2 and R 2 . As noted in [7], going to the higher dimensions or 
considering other geometries such as S d ~ 2 x R do not change the coefficients. 

But in the presence of cubic terms due to the enormous number of terms in (13.41) we need to 
consider more geometries than before. To compute the contribution of each term we consider 
cylindrical and spherical geometries inside the AdS space-time as our entangling regions (indeed 
we used 5 different geometries). Now let’s start from the following regularized metric 

ds\dz 2 + r 2 dr 2 + 1 dr 2 + (a + r n c 1 ~ n cos r) 2 ds\\ , (3.5) 

■A'Tn \ f)^ / 

where S is the entangling region (cylinder or n-sphere). 

To compute the entanglement entropy we need to find a similar relation as ED for each 
term in the Lagrangian. For example we start from the following term 

f d d x 4 /gRR 2 = n f d d xVQ7llZ 2 + 47r(l — n) f d d ~ 2 Xy/j ( 1 Z 2 + A.7 Z^n^n^) 

JMn Jm Jt. 

+ [ d d ~ 2 x^(A 1 R 4 + A 2 R 2 R 2 + A 3 RR 3 + A 4 R 2 + A 5 A' 4 ) 

J s 

+ ^ d d - 2 x^(B 1 KR 2 +B 2 RR^RP 1 '+B 3 R^K" a RP+B A RR 2 +B 5 RP a R 1 'f 3 R lia „p), (3.6) 

where, as we discussed before the second integral on the right hand side of the first line is the 
Wald contribution. To compute the unknown coefficients, ,B§, we must compute each 

value of integrals in both sides. Each integral can be written as a series, f dz c n z n . Then we 

can make a set of algebraic equations by looking at the same powers of z on both sides. You can 

see the final results of coefficients in table 1. There are some general properties in this table: 

1. The values of Wald’s entropies have different power expansion other than the value of R- 
terms or RR- terms. So the value of Wald’s entropy does not appear in our algebraic equations. 
Note that it has a non-zero value which equally appears on the integrals of the left hand side. 

2. The power expansion of A"-terms differs from AA'-terms too, but the integrals of the 
left hand side have both power expansions. Therefore we have two distinct sets of algebraic 
equations, one for Aj’s and one for B^s. 

3. There are some simple relations between the values of integrals on the right hand side at 

each dimension d due to the fact that for n-spheres the Ricci curvature is proportional to the 
metric. For example in AdSd, R ~ d — 3 and 1Z ~ d{d — 1) therefore 1ZR 2 = dRR^ v RP v or 

KR 2 = dR^n^RP and 1ZR 2 - KR 2 = d{d - 1)AT^ R^K^p. 

4. With the same reason as mentioned above the value of RR 3 and A'| are equal independent 
of dimension for cylinder or n-spheres, so we consider just one of them. It is important to 
remember that this may not be correct for more general entangling regions other than the 
cylinder or spheres. 
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Finally by solving the algebraic equations we find 


/ d d Xy/gRR 2 = n j d d x\/GJZJZ 2 
■JMn Jm 

+ 4tt( 1 - n) f d d -' 2 x^(n 2 + rni^nfn" - -K* + ^A' 2 /F 2 - ^W 2 ) , (3.7) 

Js \ 4 4 1 J 

Performing the same computations for R? and R 3 terms in the Lagrangian (see appendix A for 
a list of integrals) we find the final formula for entanglement entropy of ENMG as 

S EN MG = 2^ + AXL 2^ n yu _ 1^2 _ 3^ + _ 3 Ff 2 A 2 

+6^k 2 - i6KK lu ,w u/ + ^-n 2 - - 12^2 + im lia 'Ji„ a nfn l ')\. ( 3 . 8 ) 

Now we can find the minimal surface similar to the NMG case. By extremizing, we will find two 
solutions, one with zero and the other with a non-zero extrinsic curvature 

fi(z) = \jzl~ z 2 \ f 2 (z ) = yjzl + 2 z 0 qz - z 2 ; q 2 = -— —^_ 3/ ^°° . (3.9) 

H'Joo ^ 

By knowing these solutions we can compute the universal terms in entanglement entropy 

ci _ 27 rL(fj J f^ 0 + 2A/oo + 1) 

" 3 "" IwjZ ’ 

^ = ^y/2(l-»f* 0 )(\ + tlf 00 ). (3.10) 

If we use the relation 1 — /qq + f^X + f^p = 0 then we will observe that to have a real function 
for f 2 (z) and for positive coupling A, we must have | > A > 0 and 1+V ^~~ 3A > f 00 > 

Note that in this interval, c 2 is real valued if /oo > 1 which is satisfied automatically. We also 
have c 2 — c 2 = 3 ^ fa L (A— 2+ 3) 2 so in the region allowed by reality condition, ci > c 2 

and / 2 (z) gives the minimal entropy. For A > | the first solution is the correct one where by 
demanding a unitary dual CFT we restrict to ^ 1 " l i A ~ 1 . The different domains of validity 

of the solutions has been shown in figure Q2J 

4 General curvature square terms 

In previous section we considered a special known action of gravity and its extension in three 
dimensional space-time and we used a line as our entangling surface. To generalize the study of 
holographic entanglement entropy to higher dimensions and for more general entangling surfaces 
such as cylindrical or spherical regions, in this section we consider the general curvature square 
terms in five dimensions. We start from the following action and use a similar notation as [9] 

1 r / 12 XL 2 \ 

S = -^3 J d b Xyfg^R + £2 -I —— (XiRftvap + A 2 A^ + X 3 R 2 )J . (4.1) 

Let’s first consider S' 3 as an entangling surface, then the non-regularized bulk metric will be 

ds 2 = (dT 2 + dz 2 + dp 2 + p 2 (dd 2 + sin 2 Odcf) 2 )^ , 

8 


oM - £i i n ££ 
bEE ~ 3 e 

0(2) - ^ i n 

*EE - 3 111 e 


(4.2) 















3 


A 2- 



0 1 


II 



Figure 3: To have c\ > 0, the region I is forbidden for f\(z). This solution is valid in both regions 
(II) and (HI). The reality condition of C 2 restricts f- 2 (z ) to the region (III). In the overlap region 
III we always have ci > C 2 - 


where L = L/\J and here we have 1 — /oo + gA/^Ai + 2 A 2 + IOA 3 ) = 0 as a constraint on 
the parameters space. Doing the same steps as what we did for previous cases, one finds a the 
differential equation which gives the extremal surface for S 3 entangling surface. Again there are 
two solutions without/with extrinsic curvature 




+ 2 z 0 qz - z 2 ; 


4 - (I 6 A 1 + 22A 2 + 80A 3 )A/ oo 
— 4 + 5(4Ai + 5 A 2 + 16A3)A/oo 


Then the leading terms in the entanglement entropy will be 
:(i) _ 


qK ±j A ( Z 0 I Zq . 7T 2 L 3 (1 — 2A/oo (Al + 2 A 2 + 10 A 3 )) 

See = 4a n^2 ~ ln “ 4 -) > a i =-5- 


m'o 


(4.4) 


See ~ 4 a 2 (- 


1 


+ 


2 q 2:0 


zo 


7r 2 L 3 A§ (4Ai + 3A 2 ) 5 


y q 2 + 1 e 2 q 2 + 1 e 
For cylindrical geometry with length H we consider 


— In — + •••), 02 = — -- 

e 2 ip aJ— 4 + 5 (4Ai + 5 A 2 + I 6 A 3 ) A f 0 


ds" = 


L 2 


^ dr 2 + dz 2 + dp 2 + p 2 d0 2 + du 2 ^j . 


(4.5) 


The possible solutions are two series solutions as 
~2 


fi{z) = Z 0 - — + 
4z 0 


f 2 (z) = z 0 +qz- (1 + q )-- 1 - 

4 z 0 


(4.6) 
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with the same value for q as in equation (14.31) . Using the above solutions the leading terms of 
entanglement entropy are 


cU) _ H (A z o i z o , 
See ~ 2zo (4fll e 2 Cl n e + 


H 


7T 2 L 3 

), ci =- 3 -(l + 2A/oo(Ai — 2 X 2 — 10A 3 )), 


o( 2 ) _ 

°EE ~ 


H 7t 2 L 3 A^ 


_£0 , / £0 / 1 £0 _\ 

2,0 ilfoo (C °e 2+Cl e +C2ln e + 

Cq = —2 \J 4Ai + 3A2\/—4 + 5(4Ai + 5 A 2 + 16A3)A/oo , 
d x = —2\J 4Ai + 3A 2 \/4 - (I 6 A 1 + 22A 2 + 80A 3 )A/ oo , 

_ 8((19Aj + ( 3 IA 2 + 80A3)Ai + A 2 (A 2 + ~^A3))A/oo — f A 2 — 4Ai) 
C2 4- 4Ai J-A + 314A, -I- 3Ao + IfihUL 


• The first solutions fi(z): As we see from the results in (14.41) and (|4.7I) the leading terms 
of both entropies have the same coefficient ai independent of the topology of the entangling 
surfaces. Moreover the coefficients of the logarithmic terms (ai and ci) exactly coincide with 
the known results of central charges of the dual CFTs. In this case and in the three parameter 
family of the solutions which specify by Ag 2 , 3 , we need to have a positive value for a x , which 
restrict us to 1 > 2A/oo(Ai + 2 A 2 + 10 A 3 ). 

• The second solutions f 2 (z): In order to have a real solution in this case we must demand 
a real value for q which means that q 2 > 0 in equation (14.31) . Moreover, to have a positive value 
for leading term of entropy, for spherical region we need to have a real and positive value of 02 - 
In cylindrical case Cq must be real and negative. 


4.1 Gauss-Bonnet gravity 


To compare the domain of validity of our two different solutions, it will be easier to choose a 
specific point in the three parameter family of the solutions. The Gauss-Bonnet gravity is one 
of the interesting cases (note that it is topological in four dimensions). 

So let’s restrict ourselves to Ai = TA 2 = —4 and A 3 = 1. The constraint explained in the 
previous section then simplifies to A = —. For non-zero extrinsic curvature solution and for 


spherical entangling region, q 



Two central charges in (14.41) simplify to 


a 1 


tt 2 L 3 

^/i 


(6 - 5/oo), 


02 


4y / 27r 2 L 3 (/ 00 - 1)2 

^/oo 3 


(4.8) 


The only consistent solution for f 2 (z) happens when q is real-valued. This restricts us to foo > 2. 
In this region the leading term of S EE is positive but for S EE it is negative. The only acceptable 
solution for f] (z) is in the region where | > foo > 0 . 

In cylindrical geometry from ([4.7D and relation 1 — foo + A= 0 we have 


- 7r2jC ' 3 c f \ 

ai — 5 5 /oo), 


ci = 


e 2 L 3 2 - U 




£3 


ft 


5 


2 foo 2 

q = 


fa 


, _ — 8 V 2 tt 2 L 3 ^Jo ^1 , -^tt 2 L 3 V^ 1 (/oo- 2 ) 

c o — 7772 1 c i — « c o j c 2 — 


p?j f 2 

^ pJ OO 


/3 /*3 
^pJ oo 


(4.9) 


( 2 ) 

The second solution exists again only for foo > 2 where S EE has a positive leading order 

(1) 

coefficient but in this region S EE is negative. For | > foo > 0 only the first solution is allowed. 
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Note that in this region the central charge ci has a positive value 0. 


Cj>0 


c 2 >0 


0 


foo 


I 

3 


4 


Figure 4: The domains of validity for different extremal surfaces in Gauss-Bonnet gravity are 
two distinct regions. The blue domain for f\ (z) is restricted to | > > 0. The red region for 

f 2 (z) is bounded from below, > 2. 

As we observe from the computations in special case of Gauss-Bonnet gravity, the domain of 
validity of two solutions is independent of entangling surfaces and specifically it has two distinct 
regions. Therefore if one demands to have a unitary CFT dual theory, which can be achieved 
by restricting to the blue region in figure (J3J) then we can ignore the second solution. 

5 Quasi-Topological Gravity 

In three dimensions we studied the extended NMG as a theory with cubic curvature terms. But 
in this case since the entangling surface was a line, the value of extrinsic curvature was zero when 
we considered the first solution fi(z). Note that we can always use the first solution as a check 
of our calculations because the coefficient of the universal term is a known value of the central 
charge for the dual CFT. To study a theory with cubic curvature terms and with non-trivial 
extrinsic curvature we consider the quasi-topological gravity in five dimensions. We will show 
that while the extrinsic curvature is zero for spherical entangling surfaces it is non-trivial for 
cylindrical case, therefore quasi-topological gravity can provide us an opportunity to check the 
results of regularized integrals in appendix A. The bulk action is given by [21J 



C 2 = R, I:V aaR ,waf) - 4 R^RT + R 2 , 

£3 = R iLVp(T Rva.o(iR ( \% + \r^R^R - ^R lwpa R ,wp pR ai3 + y R^pR^R^ 

+y R naR a ^Rp M - ^R^R^R + • (5.1) 

The entanglement entropy for this case also has been studied in [9] as explained in the section 2, 
and in m using the formula found in |10) for quadratic curvature theories. Both results agree 
with the universal terms expected for this theory. 

Since we have cubic correction terms we need to use the regularized integrals given in ap¬ 
pendix A. Considering each term in Lagrangian and replacing its corresponding integrals and 
after simplification we will find the following relation for holographic entanglement entropy of 

1 For holographic c-theorems specially in odd dimensions see [22). 
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quasi-topological gravity 

S& = !/ + L 2 \(n + 7 Z» vpa n l *n p i n v j n a j - 2 ?V« + A 2 - AT 2 ) 

o r 4 

+ ^^(l5^ 2 -44^ 2 + 7^ p(T ^ 1/ ^-4(ll^ Ml/ + 37e QCTpAt ^“%-10^ pi ,^-18V^)« 

+ 8(57e pp ^ ff - 6 K^palTr - 7'Jl a l f p 'R. ua<T0 + ^KK^) («n>? - «+>(() + 2A+ + 25 A 2 
- 38A 2 A 2 + 11A 4 + 22AA 2 -104AA^A^ + 104A pa A ay A£-22AA 2 +8A^ Q A^A^ Q/3 )^ , (5.2) 

where this result includes the Wald’s entropy as well as corrections coming from the existence 
of non-trivial extrinsic curvature. Similar to GB case in five dimension we can find the value of 
holographic entanglement entropy for different entangling regions. 

• For spherical region the solutions of the corresponding differential equation for minimal 
surface are 

fi(z) = yjz o - 2 2 , / 2 O) = \JZq+ 2z 0 qz - z 2 , q 2 = -1 + (A ± ^A 2 + 3/r)/oo . (5.3) 

For each solution we find the following entanglement entropies 

See = ^ y 2 ( l ~ 6A /°° + 9 a/(L)(|I - In ^ + ■■ ■), (5.4) 

^pJOQ 

M2) _ 4tt 2 A 3 <? 4 + 2(1 - 3A/oo)<? 2 + 1 - 6A/oo + 9/r/ 2 1 2g z 0 z 0 , N 

^ (? 2 + 1 ) V2 V + ld g 2 + le n e j ' 

As we see the coefficient of the leading term (and logarithmic term) is the central charge cor¬ 
responding to the dual conformal held theory of the quasi-topological gravity. To study the 
leading order terms one may use the constraint /i/^ + A— foo + 1 = 0. To have a positive 
value for entropies both leading terms must be positive valued. Additionally we demand that 
the second solution is real valued. The positive value of the leading term in restricts us to 
0 < A < and therefore > yjj if we suppose a positive coupling. This corresponds 

to a region below the blue curve in figure [5j On the other hand this condition for leading term 
of S gives 0 < A < < A f or /^ > 1 (the regions (I) and (II) below the red 

00 ^ 

curve in figure [S]) . The reality condition for the second solution also restricts us to the region 
0 < A < 2 %~ 3 and therefore > § (the right hand side of the green curve in figure [5J. It is 

possible to show that in region (I), the leading term of is always greater than the leading 
term of , so the minimal surface corresponds to /^(z). 

• One can also look at the cylindrical entangling region. The series solutions can be found 
by solving the differential equation for f(z) 

/i(z) = £o-^d -, (5.5) 

f '2 ( 2 ) = z 0 + qz - (1 + g 2 )-— H-, g 2 = -1 + (A ± \A 2 + 3v)foo ■ 
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Figure 5: The regions allowed by positivity condition of leading terms in quasi-topological gravity 
are the regions (I) and (II). The reality condition of the second solution also restricts us to the 
region (I). 


The corresponding universal terms of entanglement entropies are 


o(i) - H (A z o 1 z o , \ 

See ~ 2-r 0 ^ ai e 2 Cl n e ^ 


ci = 


q(2) _ _ 
°EE ~ 


TT 2 L 3 
f3 f 3 / 2 

t'pjoo 

2it 2 L 3 H 1 


7r 


! l 3 


(1 - 2A/oo - 3/z/^), a, =-— (1 - 6A/oo + 9///^) 

fSf 6 / z 

^rJoo 


(<« + <4- + 4 in-). 


^0- 


~0^p/c» 2 (9 2 + 1)5 


c 0 — O' 1 + 2(1 — 3A/oo)(/ 2 + 1 — 6A/oo + 9/r/^ , — gcg , 

C2 = if 6 + ^(13 - 8A/oo)g 4 + i(16 - 20A/oo + 51/r/^)g 2 + ^(1 - 2A/oo - 3 pcfl) .(5.6) 


As we see, by comparing the coefficients of cylindrical region with spherical region, we will find 
exactly the same conditions as what we found for the spherical entangling surface. We have 
observed the same behavior when we considered the Gauss-Bonnet gravity in five dimensions. 
Therefore this suggests that at least for cubic corrections such as quasi-topological gravity (which 
the equations of motion are quadratic for AdS background) the domain of validity of solutions 
is independent of topology (here cylindrical and spherical surfaces). We must emphasis that the 
value of extrinsic curvature is non-zero for cylindrical case when we consider the first solution 
and this gives us an opportunity to check the regularized integrals in appendix A specially the 
correction parts including the extrinsic curvature. As we see the value of ci and oi exactly agree 
with known results for central charges found for the dual CFT of quasi-topological gravity. 


6 Conclusion and Summary 

In this paper we have used the procedure for computing the integral curvature invariants on 
manifolds with squashed cones which was introduced in [7j and we have reviewed it briefly in 
introduction. This method has been used already to compute the corrections to the Wald’s 
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entropy. The corrected terms are constructed out of extrinsic curvature of the minimal surface, 
and for the Gauss-Bonnet gravity, they coincide with results of Jacobson-Myers entropy. 

Here we apply this method to compute the holographic entanglement entropy for gravities 
with cubic curvature terms. Using our results we can compute the leading and the universal 
terms of entanglement entropy for theories such as extended NMG in three dimensions or quasi- 
topological gravity in five dimensions. 

In |7j the curvature squared gravities are studied and two different geometries with squashed 
cone has been used to compute the corrections to the Wald’s entropy. These corrections corre¬ 
spond to the existence of extremal surfaces with non-zero extrinsic curvature (11.41) . By going to 
the higher curvature theories such as the cubic curvature gravities which we have considered in 
our paper, one needs to consider corrections of the following form 

K\ K\, K 2 K\ KK 3 , /v 4 , n^K^K, n*KZK£, UK 2 , UK 2 , 

in addition to the usual Wald’s entropy terms. In this regard since we have ten independent new 
terms (for example see equation 13.61) then we need to consider at least five different geometries 
with squashed cone. The main difference between our work and |7j is the fact that we have 
considered the squashed cones inside the AdS d space-time where d = 5,6,7, 8 . Doing this, as 
we see from the table 1 in the appendix B, the first five terms will have different expansions in 
term of the radial direction z than the last five terms and therefore with just five geometries 
we will be able to compute the ten unknown coefficients. One must notice that in computing 
each integral (see appendix A) the integrals over regularized metric have a finite expansion 
series which must be reproduced by all integrals on the right hand side. We have computed 
different possible curvature cubed terms that appear in ENMG or quasi-topological gravity and 
results are summarized in appendix A. We have have also computed the integrals for d = 9 
(columns six and twelve in table 1) as a check for our computations. Consequently we believe 
that the coefficients of integrals in appendix A are universal and independent of the dimension 
of space-time, at least for cylindrical and spherical entangling surfaces. 

In studying the gravity theories including the quadratic or cubic curvature terms, we have 
observed that in addition to the usual minimal surfaces in the literature for spherical and cylin¬ 
drical entangling surfaces i.e. 

h(z) = \Jzl~z 2 , h(z) =z 0 -^- + --- , 
there are other types of minimal surfaces with non-zero extrinsic curvature (see figure |T]) 

f 2 (z) = \]zl + 2 z 0 qz - z 2 , f 2 (z) = z 0 + qz - (1 + d-• 

The value of q depends on the theory which we are studying and it restricts us to special regions 
of parameter space since f- 2 (z ) must be a real function. In addition to this restriction, we must 
take into account a positive value of the leading terms in the entanglement entropy as well, this 
usually is equivalent to the existence of a unitary CFT dual when we consider the first type of 
solutions i.e. fi(z). We have confirmed this result (value of central charges) everywhere as a 
check of our computations, see equations (12.111) . (13.101) . (14.41) . (15.41) and (15.61) . 

It must be noticed that although we have computed the entanglement entropy for both 
spherical and cylindrical entangling regions, but the domain of validity of the solutions are 
exactly the same and is independent of the topology of entangling surface. 

As we mentioned, the universal terms in the entanglement entropy will be proportional to 
the known central charges of the dual CFTs if we consider the first type of the solutions. On 
the other hand the universal terms of the second type will give rise to other values. Therefore 
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this question arises, which solution corresponds to the correct entropy? or if the first type of 
the solutions give the correct entropy then how one can get ride of the second type solutions? 

According to the prescription given in [1] and [2], the correct entanglement entropy corre¬ 
sponds to an extremal surface which gives the minimum entropy. To understand which solu¬ 
tion has the minimal entropy we have compared the leading order terms in the entanglement 
entropies. To summarize our results, let us compare our answers according to the order of 
curvature terms in the Lagrangians: 

• Curvature squared Lagrangians: 

By looking at the results of NMG (figure [2]) and GB (figure H|) we see that the allowed 
regions of the solutions in the one dimensional parameter space of these theories (imposed by 
the positivity of the leading terms and reality condition) are two distinct regions. So if we 
restrict ourselves to the unitary dual CFTs we can ignore the second type of the solutions. For 
these theories, this choice of extremal surfaces, has been imposed by causality constraints in 

m, 

• Curvature cubed Lagrangians: 

The results of ENMG and quasi-topological gravity have been summarized in figures (J3]) 
and ©• As we observe, unlike the curvature square Lagrangians, there are regions in the two 
dimensional parameter space where both solutions can co-exist. Interestingly, in both theories 
regardless of entangling surfaces, in this overlap region always the leading term for S ^ is greater 
than the leading term of S^ and therefore the minimal entropy corresponds to f 2 (z). So in cubic 
curvature theories we can not ignore the second solution by imposing the unitary constraints. 

Consequently, although unitary constraint (unitary dual CFT) may prevent us to deal with 
the second type of extremal surfaces ,f r 2 (z) in curvature squared theories of gravity, but it can 
not be helpful for cubic curvature Lagrangians. It would be interesting to find a constraint from 
the CFT dual side in order to restrict us to the first type of solutions fi(z), which gives the 
actual value of entanglement entropy. 
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A The curvature cubed integrals 

In this appendix we have summarized the value of curvature cubed terms which appear in 
different Lagrangians 


/ d d Xy/gRR 2 =n d d x\JQTZTZ 2 
Jm„ Jm 

+ 4vr(l -n) f d d ~ 2 x^(n 2 + - ^/v 4 + -K 2 K 2 - , (A.l) 

J S — 



= n d d x\[QR? + 127r(l 

Jm 


n ) J d d 2 x - v /7 1Z 2 , 


(A.2) 



d d Xy/gR 3 


= n d d xy/GlZz + 67r(l 

Jm 


n) [ d d 2 x y /^{TZ p a n ua n p n u i 
J E 


KK pv W 11 ') , 

(A.3) 


[ d d x^fr lupC7 R„aopR a f P = n[ d d xVQn^ pa n ua<yp n a ^ p 
JM n Jm 

+ 67t(1 - n) f d d ~ 2 xy/j (R afl,3f> 'Ti P upa [n”n rTl n ol jn j p - n\ripn a jn aj ] - ^K 2 + ]-K± 

J s 

+ - yRK 2 - , (A.4) 


[ d d x^RR pupa R pvp ° = n f d d x 4 fgRR, IIL , prJ R piJp(J 
JMn Jm 

+ 4vr(l -n)J d d - 2 Xy /7 ( [R, wprj R pvpa + 2RR^ prr n pi rj p n uj ni] - K 2 K 2 + K 2 - 2AA 2 ) , 

(A.5) 


[ d d Xy/gRa R P v P (jR ilvpa = n [ d d x^R 0 °R pvprJ R^ pa 
JMn Jm 

+ 27t(1 - n) [ d d ~ 2 x^y(2R pi/pa Ra [n pi n p n uj n J a - n^n^n^nU + R pup(J R pvpa n a in 0 
Jt, v 

- K 2 + K 4 - ^KK pi/ TZ lw - ^ Kpa n™K p + | RK 2 - ^K pa K^n pua p 


(A.6) 


/ d d x 4 jJ]R va R pp R pi y P a = n / d d x^R uu R pp R pupa 
lM n Jm 


+ 27t( 1 - n) J d d 2 Xyf^(2R p , vpa R pp n v i n ai + R pp R ua [n fli n l p n vj n : l - n pi vJ a n vj n j p \ 

23 7 12 

- K 2 K 2 + A' 2 - —KK pv n pv + —K pa R av K p - —KK 2 + -K pa K^K pua p 


(A.7) 
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B Table of Integrals 


J d d - 2 xj 7 --- 

AdSg 

AdSf 

AdSf 

AdSf 

AdSf 

AdSf 

AdSg 

AdSf 

AdSf 

AdSf 

AdSf 

AdSf 

A 4 

2 

64 

162 

2U48 

3 

625 

6912 

5 

0 

0 

0 

0 

0 

0 

k 2 i< 2 

2 

32 

54 

512 

3 

125 

1152 

5 

0 

0 

0 

0 

0 

0 

kk 3 = K‘i 

2 

16 

18 

128 

3 

25 

T9? 

5 

0 

0 

0 

0 

0 

0 

K, 

2 

8 

6 

‘S2 

3 

5 

o2 

5 

0 

0 

0 

0 

0 

0 

77 K 2 

0 

0 

0 

0 

0 

0 

-40 

-320 

-540 

-1792 

-1400 

13824 

5 

KK I 1 U W‘ AJ 

0 

0 

0 

0 

0 

0 

-8 

-64 

-90 

-256 

-175 

_158b 


0 

0 

0 

0 

0 

0 

-8 

-32 

-30 

-64 

-35 

25b 

5 

77 K 2 

0 

0 

0 

0 

0 

0 

-40 

-160 

-180 

-448 

-280 

23U4 

5 


0 

0 

0 

0 

0 

0 

0 

-8 

-12 

-32 

-20 

-32 

I d d Xyfg--- 
JM 

AdSg 

AdSf 

AdSf 

AdSf 

AdSf 

AdSf 

AdSg 

AdSf 

AdSf 

AdSf 

AdSf 

AdSf 

RR-2 

0 

32 

108 

512 

500 

1152 

-80 

-640 

-1080 

-3584 

-2800 

27648 

5 

R 6 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Rs 

0 

0 

0 

0 

0 

0 

-48 

-384 

-540 

-1536 

-1050 

921b 

5 

R^ pa R„ a apR a f P 

0 

24 

36 

96 

60 

96 

12 

0 

-36 

-128 

-90 

768 

RR^ P aR pvpa 

0 

64 

144 

512 

400 

768 

-320 

-1280 

-1440 

-3584 

-2240 

18432 

5 

R Wvpoc d a 

XL p,ispcr- LL 

0 

16 

24 

64 

40 

64 

-80 

-384 

-384 

25t)U 

3 

-480 

_ 3584 

5 

/? R^ P R VG 

J - L p,upa J - L 

0 

32 

72 

256 

200 

384 

-28 

-224 

-324 

281b 

3 

-650 

-1152 

X 

7THz.fl' 2 

rl/2 

nzjo6 

7T 2 

a 

7T 2 L 

7 v 3 aL 2 

7 r 3 a 2 L 3 


*1/2 

Tlfoi 

7T 2 a 

7 r 2 a 2 L 

7 r 3 a 3 L 2 

7 r 3 a 4 L 3 

Za 3 

La 2 


Z^foc 

q ,3/2 

Laz 

Lz 

z 1 


Z 4 foo 

= *3/2 

Z 5 fcJ, 


Table 1: The value of each integral can be found by multiplying its value to a proper coefficient 
shown on the row by x sign. For example: fg 2 (Px^fK 4 = QAn z ff 2 / La 2 . For the second part 
of the table each value additionally must be multiply by tt( n — 1). 
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